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TECHNICAL MEMORANDUM 


AN EXPERIMENTAL INVESTIGATION TO DETERMINE INTERACTION 

BETWEEN ROTATING BODIES 

(MSFC Center Director’s Discretionary Fund Final Report, Project No. 279-00-16) 


1. INTRODUCTION 


Present physical theories, both phenomenological and quantum, while in a broad agreement 
with experimental evidence, on a deeper level of logical analysis, encounter many inconsistencies or 
problems. Although these problems have no immediate practical consequences, they stimulate the search 
for new fundamental models of reality. Very active areas are now superstring theories, geometrodynamic 
extensions of general relativity, such as Kaluza-Klein theories, topological geometrodynamics theories, 
and twistor or torsion field theories. This new theoretical breeze brings new challenges to experimentalists. 

As our quest for better understanding of the physics governing nature continues and more 
experimental data became available, we begin to see some remarkable similarities between the macro- 
cosm and microcosm. Perhaps the most important is the rotational motion of matter at any length scale. 
Planets, stars, galaxies, and clusters of galaxies are all spinning objects, and at the microscopic level, all 
the elementary particles (=1,000) are spinning entities. Consequently, we witness that matter spins at all 
scale levels. If one assumes a spinning distribution of electric charges than are in accord with the law of 
induction, a magnetic dipole will result. Although such a classical model cannot be directly translated 
into the microworld, the general cause for the phenomenon of the magnetic moment remains the same. 

The magnetic moment of the particles is further related to its spin through the gyromagnetic 
constant, and theoretical predictions for this parameter are in excellent agreement with experimental 
data. However, there is mounting evidence on both the microlevel and macrolevel that our physical 
model of spinning objects needs to be reassessed. Many of the theories mentioned above and models 
of particles as nonsingular extended objects, usually based on classical electrodynamics, point to the 
existence of some new types of interactions. Some of these theories, while maintaining the essence of 
the classical physics as a subset, predict interaction between spinning macroscopic bodies. Also, reports 
exist of experiments, conducted over the last hundred years, that point to such phenomena. It is difficult 
to assess the credibility of these data since others have not repeated the experiments. 

The new paradigm in physics revolves around the concept of vacuum and/or multidimensional 
geometry. Our commonly accepted understanding of the vacuum is based on the Dirac sea of particles. 
Here the energy spectrum of the particles has a gap and the lower energy band is almost fully occupied. 
When a particle acquires sufficient energy, it is transferred from the lower to the upper energy band. A 
free particle is created this way, and a hole is left in the lower band that represents an antiparticle. From 
the quantum principle of uncertainty, there will always be some fluctuations in the system. This means 



that the vacuum actually constitutes a complex fluctuating field of various particles. Wheeler 1 evokes 
that this particle flux has an energy density on the order of 10 94 g/cm 3 . The Casimir force, recently 
measured, 2 that acts between two electrically conducting bodies confirms the viability of a vacuum 
as a real medium. In short, the physical properties of a vacuum depend on the presence and strength 
of macroscopic electrical, magnetic, and gravitational fields. 

Spin is the fundamental property of elementary particles. For example, neutrons, electrically 
neutral particles, have a spin of til 2, as do electrons. Many of the modem quantum and classical theories 
also consider the physical vacuum as having spin or torsion properties. If so, then some macroscopic 
manifestation similar to the Casimir force would be present. Modification of the local vacuum spin 
properties by electromagnetic fields or by rotating bodies is then possible. However, no a priori 
estimates of the effect can presently be made because the current vacuum models are inadequate. 

Scientists in the 19th century liked to consider the free space between bodies as a medium with 
fluid-like properties so spin can be transmitted. Beginning early this century many investigations have 
been conducted with the intent of establishing an influence of solid body rotation on electrical, mag- 
netic, and gravitational properties. Einstein, recognizing that rotational motion was not correctly treated 
in his first general relativity theory, proposed, with Cartan, the first gravitational theory of spinning 
bodies in 1 923. 1 Briefly, this general relativity theory, in its classical formulation, predicts a small 
correction to the Newtonian gravitational force for spinning macroscopic masses. This correction is 
rather difficult to detect under presently achievable rotational speeds. Several new proposed superstring 
models dealing with reconciling gravity and quantum theory point to bigger effects. Also, twistor or 
spinor theories, comprising natural extensions of the general relativity, assume the existence of 
additional forces between the spinning bodies of unspecified magnitude that can only be determined 
experimentally. Some experiments were reported decades ago with no definite claims, and only a few 
measurements have been recently attempted. Also, we want to emphasize that there are several ways 
rotational motion can be included into the general relativity theory, as attested to by many later models, 
that address the necessary existence of gravitomagnetic fields. 3-8 

Although clearly overshadowed by theory, there have been some limited experimental works 
conducted in this area. Myshkin 9,10 conducted a series of experiments purported to discover a relation- 
ship between macroscopic fields and rotation. Experiments utilizing gyroscopes that have not received 
much exposure were later performed in the U.S.S.R. by Kozyrev. 1 1 Other obscure Russian literature on 
experiments and theories related to non-newtonian mechanics includes the following: 

• Spartak Poliakov’s book with experiments on complex gyroscopic systems that exhibit 
non-newtonian behavior. 

• A theoretical book by Gerlovin, whose approach is based on Einstein’s general relativity 
theory, that extends the dimensions of the world from four to seven. He obtains a spectrum 

of the elementary particles, rest masses, and life times that agree with existing data reasonably 
well (six digits or so). His approach is partially based on his discoveries that some classical 
orbits of the point-like charges could be nonradiative; similarities with the string theories are 
also evident. He also predicts some unknown gravity-related phenomena. 
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• Veinik’s book addresses thermodynamics of real processes. He describes his experiments 

on gyroscopes as well as on torsion pendulums, the effects of which are not easily explainable 
by existing mainstream physics. 

• Shipov’s book relates his ideas on torsion fields. Mathematically his book is sound; whether 
his theoretical model is applicable to the real world remains to be seen, but his ideas attract 
some of the best scientists in the field. 

A detailed experimental study that suggests a possible force between spinning bodies has been 
conducted by two physicists from Japan and reported in great detail in Phys. Rev. Lett. 12 In short, a small 
incremental decrease in weight was detected by a chemical balance as a function of increasing the 
gyroscope’s rotation rate in a clockwise direction. Curiously, there was no effect observed when rotating 
counterclockwise. Subsequent experiments conducted in the United States, 13 utilizing a very precise 
gyroscope, did not find any weight changes. 

A French group 14 did record a change, albeit 20 times smaller than that observed by the Japanese 
group. Here it is relevant to state that the United States and French teams used different types of gyro- 
scopes and did not accurately reproduce the Japanese experimental conditions. It has been conjectured 
that the effects observed by the Japanese were due to vibrations of the gyroscope rotor that utilized 
average quality ball bearings. 15 In view of this premise, it is now interesting to compare the results with 
those of the Russians 1 1 who, independent of the Japanese, also used a chemical balance to measure 
gyroscope weight as a function of rotation. Initially, to within six significant figures, they did not see 
any weight changes. However, upon intentionally introducing small vibrations into the shaft, the weight 
of the gyroscope was affected in the fourth digit and only in one rotational direction, as in the Japanese 
study. The gyroscope was suspended on a rubber strip to suppress propagation of vibrations to the 
balance; this provision did not affect the observed phenomenon. In support of the above observations, 
Shipov 16 theoretically shows that nonsymmetric rotation is necessary to generate a change. 

The same Japanese team recently measured the free fall of gyroscopes in a vacuum tube. They 
again observed a reduction of the acceleration for clockwise-rotating gyroscopes. There has also been 
one report (apparently not refereed) of an experiment conducted in the United States 17 of freely falling 
gyroscopes in which the authors claim to observe changes in drop velocity due to the imposed rotation. 

The above studies utilized measurements based on rotating objects. Consider now the series of 
very well-planned experiments (that passed considerable critical review) elegantly conducted by Brush 
in the early 1920s. 18-20 After an introduction that evokes the essence of string theory, he considers 
simple motion of pendulums and drop velocities of different metals (zinc, bismuth, and lead) fashioned 
to have the same weight, dimensions, and centers of mass. Invariably the denser metal swung or dropped 
faster, a direct contradiction to the established laws of physics. He concludes that there is something 
intrinsic to materials themselves, making particular note of the peculiar behavior of bismuth. 

Bismuth is not a typical metal. In contrast to most, it expands upon freezing. Bismuth is the most 
diamagnetic of all metals and only mercury has a smaller thermal conductivity. It has high electrical 
resistance and the greatest increase in electrical resistance when placed in a magnetic field; i.e., magne- 
toresistance. Finally, bismuth has 209 nucleons, the largest of any stable element. These characteristics 
bring to mind Wallace’s detailed patents. 21-23 In short, Wallace, while at General Electric, invented 
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a machine that spun a mass to promote nuclear spin alignment. This, in turn, generated a secondary 
gravitational field that he called a kinnemassic field; i.e., a gravitomagnetic field. It was necessary, or 
perhaps optimal, for the spinning mass to have an odd number of nucleons. Wallace employed a brass 
wheel, although he notes bismuth to be preferable. To the best of our knowledge, outside of Wallace’s 
results, there were no subsequent studies (at least published) confirming or denying his claims. This has 
been a point of numerous speculations. 

We note that many researchers, apparently independent of each other, have conducted 
experiments that imply the existence of unusual forces due to rotation and/or electric/magnetic field 
interactions. These efforts are still met with some skepticism by the scientific community. We have, 
consequently, fashioned our initial experiments after them. 
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2. INITIAL EXPERIMENTS 


Our preliminary experimental work is patterned after some earlier efforts, primarily those 
conducted in Japan that reported a small incremental decrease in a gyroscope’s weight as a function 
of increasing its rotation rate in a clockwise direction, 12 and those of Wallace 21-23 conducted in early 
1970s that consisted of spinning a brass wheel. 

In contrast to the Japanese experiment that utilized gyroscopes, our work consisted of placing a 
small, high-speed electric motor directly onto the plate of a microbalance. The general setup can be seen 
in figure 1. Our experiments took a number of progressions. Initially, to reduce vibrations, insulating 
foam was wrapped around the motor. Shortly after initiating rotation (up to 30,000 rpm), a reduction 
of weight was detected that, upon stopping, did not return. Initially surprising, we soon came to realize 
the heat of the motor had vaporized some of the insulation with which it was in contact. Subsequently, 
the insulation was removed and the motor suspended in a plastic laboratory bottle. When the motor was 
spun in this configuration, a small decrease in weight (=20 rfig) was again detected. In contrast to our 
first attempt, the weight returned after rotation ceased, and the result was reproducible. However, upon 
sealing the motor in the bottle (fig. 1), initiation of rotation produced no decrease in weight. A simple 
calculation showed that if heating from the motor slightly increased the local air temperature, the change 
in density could account for the observed weight loss. The above experiment is not exactly like the 
Japanese effort but is valuable in pointing out the care needed in taking and, in particular, assessing 
measurements. 



Figure 1. High-speed motor assembly as it sits on the microbalance. 
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A mechanical, air-driven device to promote rapid rotation (up to 20,000 rpm) of a metal disc 
(in this case bismuth) was also constructed (fig. 2) and used to conduct experiments. The first experi- 
ment consisted of spinning the disc in the presence of the Hall probe with the intent of detecting a change 
in the adjacent magnetic field; none was detected. The second experiment focused on determining if 
rotation of the disc would influence movement of another disc. Here, a second bismuth disc (fig. 2) 
was suspended <1 mm above the rotating disc; a plastic film was placed between the two to elimi- 
nate convection effects. A mirror was attached to the suspension wire on which a laser beam was shown 
and reflected onto the adjacent wall. Thus, a very sensitive means of detecting any motion of the 
suspended disc was ensured. The lower bismuth wheel was then accelerated to high rpm while the 
reflected laser beam was observed. After numerous tries and variations, it was concluded that the lower 
rotating wheel did not influence the upper suspended wheel. The results, or perhaps more accurately 
nonresults, led to our third experimental investigation in which a spinning brass disc was utilized. 



Figure 2. Air-powered device designed to rotate a metallic disc, in this case 
bismuth, at high speeds. Another bismuth disc is seen suspended 
above the rotating disc. 
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3. THEORY OF THE NONAXISYMMETRIC FARADAY DISC 


The Homeopolar machine, or so-called Faraday disc, is perhaps the simplest example of a 
spinning body in presence of a magnetic field. It consists of a metallic disc with a static magnetic field 
applied along the axis. When rotated, a voltage across the radius of the disc is generated. This was the 
crucial experiment, performed more than 1 50 yr ago, in discovering the induction principle of electrody- 
namics. Since then many attempts were made to refine our understanding of its underlying mechanisms, 
and discussion of this matter is still continuing. 

One of the unresolved questions is whether there is an additional degree of freedom for the 
magnetic field. Consider, as an example, an axisymmetric magnetic field induced by a permanent 
magnet. Now, let the magnet spin around its axis. This prompts the question of what happens to the 
magnetic field. Do the magnetic field lines rotate as well, or is it just magnet motion that has no conse- 
quence to the magnetic field? Several clever experiments have been proposed, but doubt still percolates 
through the relevant literature. Consequently, there is growing interest in fundamental research on 
possible unknown effects induced by spinning of macroscopic bodies. For example, we note the torsion 
balance experiments by the Eot-Wash Group at the University of Washington 24 looking for new short- 
range interactions, existing in accordance to some recent theoretical speculations. Gravity probe-B 
experiments by NASA will accurately measure rapidly spinning spheres in space with the aim of 
investigating gravitomagnetic effects predicted by Einstein’s general relativity. 

Our initial interest, which led to this work, stemmed from the ideas of Wallace that were presented 
in a series of the U.S. patents. 21-23 Wallace, while at General Electric, investigated a rapidly spinning 
brass wheel and conjectured the existence of a novel field. No further work was done to confirm or deny 
his claims. Experimentally, a brass wheel 8.6 cm in diameter and 1.88 cm in thickness was rotated 
up to 28,000 rpm and a faint signal measured by a Hall probe was detected at a distant point. In an 
attempt to explain his results, he pointed to spin interactions. A convenient way to describe the direct 
spin-spin interactions between the spin-polarized beams of particles and the target is by introducing the 
concept of pseudomagnetism, an idea similar to the Weiss molecular. Our hypothesis is that pseudomagnetism 
can explain Wallace’s observations. 

In our setup, we avoided the use of any magnetic material and spun a 4-in-diameter, 1 .4-in-thick 
brass wheel by an air turbine. Furthermore, we shielded our system from interference with Earth’s 
magnetic field (fig. 3). 
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Figure 3. High-speed brass wheel assembly used in this investigation. 


To detect a possible but unknown field (gravitomagnetic), we used both the Hall probe and a 
much more sensitive giant magnetoresistive sensor. Within our sensitivity limit (10 -4 G) and a rotation 
rate of 8,000 rpm, no signal was detected. However, with no magnetic shielding, a strong signal, its 
magnitude increasing as the rotational speed increased from zero, was measured at the rim of the wheel. 
Unexpectedly, the signal went through a peak value and began dropping down as the rotation rate 
increased further, even passing through the zero value. For this effect to happen, the Earth’s magnetic 
field would have to be orthogonal to the axis of the wheel. This configuration is nontypical for Faraday 
disc studies, and no relevant literature was found. 

In an effort to understand the phenomenon within the frame of magnetohydrodynamics, a theo- 
retical model of a very long cylinder was developed that qualitatively reproduces the effect. In essence, 
the explanation is as follows. Consider an arrangement of a wheel with its axis positioned horizontally 
in a vertical magnetic field. When the wheel begins to spin, an electric current is induced in it that, in 
turn, induces a magnetic field that is in a horizontal direction. Now, the total magnetic field will be 
in a nonvertical direction. As the speed of the wheel increases, the induced magnetic field will become 
nonhorizontal. Specifically, it will have an upward vertical component. Now the induced current will be 
orthogonal to the sum of the external and induced magnetic field, triggering further rotation. This self- 
consistent treatment will result in an upward drift of the direction of the maximum of the induced 
magnetic field. To conclude, as the speed increases, the magnitude of the induced magnetic field 
increases and the direction of the maximum shifts upward. The detector placed along the horizontal line 
crossing the wheel axis will initially display a signal buildup due to an increase of the magnitude of the 
induced magnetic field. Then, when the magnitude of the field increases, the maximum shifts from the 
horizontal to the vertical position, resulting in a signal drop. 
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Here we would like to comment on the Wallace experiments. The faint signals he observed 
could result from a lack of magnetic shielding. Use of a steel base and support bars is a potential 
source for measurements to be contaminated by parasitic magnetic fields. Although we did not reach his 
28,000-rpm level, our disc is markedly bigger and the detection system is more sensitive. Moreover, 
the sensor is based on alignment of the nuclear spins, and should be preferable to the measurements 
of spin-spin interactions. 

In this work, we developed a theoretical model of a nonaxisymmetric Faraday disc. First, we 
formulate an analytical model of an infinitely long, rotating cylinder in the magnetic field orthogonal to 
its axis. This model has its own merit as it clearly demonstrates the basics of the observed phenomena. 
The second is a finite cylinder model. Essentially, it is a semianalytical model. However, it allows 
for very good approximate analytic solutions. The presented solution could, for example, be used 
as a benchmark for numerical magnetohydrodynamic modeling. We note here that the considered model, 
due to its simple geometry, is a basic magnetohydrodynamic problem related to dynamo theory. Finally, 
in the last section, we present our experimental data and compare it with the proposed theory. 
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4. LONG NONAXISYMMETRIC FARADAY DISC 


In this section we develop a simple analytical model of an infinitely long, rotating metallic 
cylinder in the presence of a uniform magnetic field orthogonal to the cylinder’s axis. 

The basic equations of this magnetohydrodynamic model are well known and are given for 
completeness in section 5. As a natural choice, we utilize a cylindrical coordinate system, with the wheel 
axis corresponding to the z axis; we do not have the z component of the induced magnetic field. Further- 
more, an angular dependence is harmonic with only one basic component. This allows one to introduce 
a complex representation of the magnetic field that could be related to the experimental value through 
this particular choice: 


fi( exp ) = Re[B(r)exp/<p] . 
For the r component, the governing equation reads 


3 2 3 3 

— t + -— 




V 


3r 


dr 


B r = /'jU<7£y(fl 0 + B r ) , 


( 1 ) 


( 2 ) 


where Bq is an external magnetic field, co is the angular velocity, // is the magnetic permeability, 
and <7 is the electrical conductivity. 

The solution of this equation can be expressed through Bessel functions and is 

B r =-B 0 +CJ { {qr)/r , (3) 

2 

where q = - i/ioco , and the corresponding expression for the B ' component is 

B (p = ~iB 0 - iCJ 1 (<7'‘) / r + iCqJo(qr) . ( 4 ) 

In the free space around the cylinder, the induced field is given by 

5 space = c7r 2 R space = _ jC f ; .2 (5) 

In this work, we do not consider magnetic materials, so the magnetic permeability is assumed 
uniform throughout the region. This implies that the magnetic field is continuous across the boundaries 
and we obtain for unknown coefficients C and C 
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( 6 ) 


C = 


2Bp 

qM<? r o) 


C' = -Bo r O 


J\M 

qr Q Jo{qr 0 ) 


Then for the value of the r component of the induced magnetic field at the surface, we obtain 


fi r( r 0) = - fi 0 + 


2BMqr 0 ) 

<?^W r 0 ’ 


(7) 


Utilizing equations (7) and ( 1 ) we can plot the induced r component of the normalized magnetic 
field, at the surface as a function of the angular coordinate for a set of q values. The corresponding plot 
is presented in figure 4. 



Figure 4. Angular dependence of the r component of the normalized 
magnetic field. The field is calculated at the surface of the 
disc for a set of 11 q values. 


A plot of the r component of the magnetic field at various angular positions as a function of 
angular velocity is depicted in figure 5. This graph can be directly compared with the experimental data. 
The bottom curve corresponds to the direction of the external magnetic field, cp= 0. Then the sequential 
values for the angles are 0.5, 1, 1.5, 2, and n (in radians). The value of 1 .5 corresponds to the horizontal 
sensor position and was used in most of the experiments. 


11 




1 



q (radians) 

Figure 5. Dependence of the r component of the normalized magnetic field 
on the skin depth parameter, q, for several angular positions: 

0, 0.5, 1 , 1.5, 2, and p, in radians. 

The presented model should give practically correct values at the center of the cylinder rim of the 
induced fields for the case of aspect ratio on the order of 10 or larger. Physically, for the high rotational 
rates, the fields at the center of the rim should also be given correctly by this model as the skin depth for 
the fields becomes small. Thus, the fields at the sidewalls will not overlap and will not penetrate to the 
middle zone. In other words, sidewall effects will not be present at this position. Plots of the magnetic 
field lines for different values of the skin depth parameter, q, are presented in figure 6. For small rota- 
tional speeds, the magnetic field penetrates through the metallic cylinder. As the speed increases, the 
lines are repelled from the center of the cylinder. Also, rotation of the lines is evident. 




Figure 6. The streamlines of the total magnetic field for a set of the q values: 

1, 3, 5, and 10. The external magnetic field goes from left to the right, 
and the wheel rotates counterclockwise. As the speed increases, 
the magnetic field is repelled from the wheel. 
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5. FINITE CYLINDER MODEL 


In this section, we develop a theoretical magnetohydrodynamic model of a nonaxisymmetric 
Faraday disc of a finite length. We consider an electrically conducting solid cylinder with electrical 
conductivity, o; radius, r 0 ; and height, 2 L. The cylinder rotates around its axis with the rotational 
velocity, co. An external uniform magnetic field of strength, B 0 , is applied in a direction perpendicular 
to the cylinder axis. The geometry of the model as well as the cylindrical coordinate system used here 
is depicted in figure 7. We further assume that the external magnetic field is turned on sufficiently slow 
and the mechanical rotational accelerations are small. For this stationary case, we have 

VxB = juJ, J = cx(E + vxB), V xE = 0 , VJ = 0, VB = 0, ( 8 ) 

where B is the magnetic field, J is the electric current density, v is the velocity of the medium, and E 
is the electric field. 






Figure 7. Schematic of the geometry and the cylindrical coordinates. 


13 



First, we eliminate the current from the equation for the magnetic field and obtain 


V x B = in t(E + v x B) . 

Further, in our case, the electric field can be expressed through the scalar potential, O, as 


E = -Vcp. 


( 9 ) 

( 10 ) 


The required velocity of the rigidly rotating cylinder is 

v = care^ . ( 1 1 ) 

At this point we introduce, for greater convenience, complex fields according to the following rule 
between the experimental B (ex P) and complex B fields: 


B( ex P) = Re[B(r)exp/<p] . 


( 12 ) 


Then the externally applied uniform magnetic field, in its complex representation, is 

B( ext ) _ # 0 e 7 . + iB^ty , 


(13) 


and we have 


vxB (ext) =-tyr5 0 e, . ( 14 > 

Since the problem under consideration is linear, higher angular harmonics are not involved, and 
thus the angular derivatives can be replaced by the imaginary unity number JL _ ; 

d<p 

Also note, we are solving for the induced magnetic field that evidently drops to zero at distances 
sufficiently far from the spinning wheel. We scale distances by the radius of the cylinder and the electric 
potential by (Or 0 2 . 

After simple algebra, we obtain from equations (8)-(10) the following set of governing 
equations: 


1 d 1 

+ -- ^ + 


dr 2 ' rdr ,- 2 dz 


B-+q 2 B,= 0 


(15a) 
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where q 2 = -ijiooor 2 , 


dr 


1 3 1 

+ ~—-^r + - 


dr 


dz* 


® + q z ®=2B z 


dB r dB, q 2 dBfp / . 2 (d<$> 

dz dr r dz r V. dr 

7 f- r ( rB V )~~ B r = ~ ic l 2 {^ + rB 0 +rB r) ■ 


V dz 


A genera] solution for equation (15a) is a superposition of the following functions: 


(15b) 

(15c) 

(15d) 


B z = C\ exp(ez)/j (£r) , k = -Je 2 - ifioa) = ^|e 2 ^q 2 , ( 

where C j and e are nonspecified at this moment. 

This form, in turn, leads to the following solution for the potential, O, and the rest of the field 
components: 


(Or 


d> = CV\(Ar)-C, — J'Akr) 
k 


(17a) 


/ 

B r = C, —J.(kr) + C. — J'Akr) 
‘ ecor k 


(17b) 


B v = c i^JA ki A +c -— J 'A kr ) 

v k'r £(0 


(17c) 


Now, in order to obtain the required coefficients, we need to consider boundary conditions at the 
surface of the wheel. As the surrounding cylinder medium (air) is an electric insulator, the normal to the 
surface component of the electric current has to vanish. At the rim (/- = 1 ), we should then have 


~+rB,=0, 

dr 


which gives us the relationship between C x and C 2 as 


C l coJ l {k) = C 2 k*J' l (k ) 


(18) 


(19) 
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In principle, we have three sets of possible solutions: 


( 1 ) 


Ci=C 2 


k? /[(*) 

0) J\(k) ' 


(2) C t =0, J\{k) = 0 . 

(3) C 2 = 0, J,(k) = 0 . 


( 20 ) 


Cases (2) and (3) yield B,(r= 1)=0. 

At the sidewalls (z=±L), we have from the condition of/, being zero there: 


30 

3z 


= -rB 0 


-rB,. 


( 21 ) 


Finally, in order to write down the components of the magnetic fields as series expansions, 
we note that as a function of z, B,., and By are even, and O and B. are odd. 

Below, we write down our general solution in terms of these three series expansions as outlined 

above. 


First consider case (1 ). As the z component of the magnetic field is an odd function of z, then 
it is convenient to expand it into a Fourier sine series as 


B 




MV) 

(An) 


sin p n : 


(22a) 


, -,2 2 2 7 T ( 

where A n = p n -q , Pn=~\ n 

The rest of the fields are 



n = 1,2,3... Also note that cos(p n L) = 0 and sin (p, ,£) = (-!) 


/i+l 


Bi 2) =0 


(22b) 




Pn 7 |(V) , ? 2 / l'(V) 

/ \ I . / . 


M r h{K) XnP,A{K)\ 


cos p„: 



-E*>, 


<? 2 / l(V) + Pn jtV) 
P/i A 5 / i / ( A ») r 7 i( A j 


COS /? /f Z 


(22c) 


(22d) 


16 



* ( 1 ) = 5 >„ 


/ h{K r ) ! r /f(V) ^ 
A « 7 l( A «) 


sin(p„z) 


(22e) 


Consider now case (2), corresponding to the solution with no z component of the magnetic field. 
This set also includes the z independent term. We have 


flp) = 0 

(23a) 

»(2) _y. iq 2k h J\{ K n r ) cos h(e„z) 

9 ' l £ n J \{ K n) cosh ( £ n L ) ^ 

(23b) 

nt2l_Y * J \{ K n r ) cos h( e „ z ) „ 

' - lA » £ „ f y lh )cosh( £ „t) 

(23c) 

A(w) sinh(e„z) 
n J\(k„) cosh(e /? L) 

(23d) 


where we have 


Bor ~ ~B() + CqJ l ) / r , B()(p ~ -‘Bo + *Q) qJ{ ( qr ) 

From the boundary condition equation (21) at z=L, 


K n J\( K n r ) 


£,i M K n) 

This equation can be readily solved yielding coefficients 


+ C 0 J l (qr) = 0 


A„ = ~- 


2C 0 E„qJ{ (q)Ji (k„) 


Note that B 2(p {r,L) = B 2(p {\,z) = B 0(p {\) and B 2r (r,L) = iB 0(p {\) . 


(24) 


(25) 


(26) 
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Consider now a set corresponding to the case (3). Let 

J i (P n ) = 0 , Yn=i A?-^ 2 • 
Then the associated solution of equation (15) is 

fi (3) = y c MM sinh(y„z) 

2 7 " MPn) cosh (y n L) 

o(3) = y c '7/; cosh(y /? z) 

9 „ " P 2 n r MPn) ™sh(y n L) 

fi (3) = y c Yn MM cosh(y /t 2 ) 
r 7 n fi n J 0 ((3 n ) cosh(y /( L) 


( 27 ) 


(28a) 


(28b) 


(28c) 


<2,(3) _ _y c r 7i(^„r) sinh(y /t z) 
^ " A, cosh(r„L) ' 


(28d) 


At this point the boundary conditions for the electric current are fulfilled exactly, but the solution is still 
not fully specified, as the coefficients C n and D n are arbitrary and not defined. The remaining condition 
on the solution of our problem comes from the far-field condition of the vanishing solution at infinity. 
This condition is nontrivial, as it requires formulating a solution in the free region outside the medium 
and matching it at the boundaries. In order to avoid this, a surface integral formalism can be invoked 
instead. Indeed, outside the medium, the magnetostatic equations can be used and thus the theory 
of a scalar potential can be applied. We consider a potential, u, that satisfies the Laplace equation: 


u(r)= u(p,z)e l<p , V 2 « = 0 . 

Then, in free space, the magnetic field components can be represented as 


(29) 


du 


-.'I-. 


“ = -ipB<p • 


From the Green theorem. 


^«(PO’-0> = ^ 


3u(p,z) 

dn 



— A-*»)rfs . 

r-r 0 


where Q=4n if / 0 belongs to the region surrounded by the surface, 5, otherwise Q is zero. 


(30) 


( 31 ) 
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Integration over the angular variable can be readily performed, leading to the cylindrical Green 
function, which can be represented in this particular form: 


G = 27r\j l (kp)j l (kp 0 )e k \ z Z °U . (32) 

0 

For our geometry, we obtain from equation (3 1 ) 


^“(po z o)=- j dzr 0 


-L 


±u(pz)-u(pz)^ 


(p z PO z o) p=r 0 


>0 


0 3 ! 

+ jpdp —u(pz)-u(pz)— G (pzp<) z o)|z=-L 


3z 


dz 


r ° r 3 3i 

-\pdp\ ^u(pz)-u(pz)— G(pzp 0 z 0 ) \ Z=L . 

r\ L - 


(33) 


In terms of the magnetic field components, we have 


B,.(\z) G (lzp 0 z 0 ) + iB 9 (lz) ~^{P=^ z Po z o) 

1 

+ \pdpB z (p, -L)[G(p, -L, p 0 z 0 ) + G(pLp 0 z 0 )] 

0 

+ ij p 1 dpB (p (p, - L ) j- [G (p, z = -L,PqZq)-G (p, z = L, p 0 z 0 )] , (34) 
0 z 

where the parity properties of the fields were already taken into account. In the above expression, 
the observation point r 0 is outside the medium region. When it approaches the surface from outside, 
the condition for the unknown coefficients C n and D n results. Below, the two surfaces of the cylinder 
will be treated separately. We first consider the sidewall. In this case, z 0 =-L-e, 0<r 0 < 1 , where £ is infini- 
tesimally small. Using the Bessel function representation for the Green function, equation (32), we 
can put equation (34) into the following convenient form: 


-4niPo B (p(pO z o)=-j dz 
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CO 


-ipoB^poL) = J dkJ x (kp 0 )U{kL) , 


(35) 


where 


U{kL) = - J dz\B r {\z)J x {k)+ iB^lz) kJ{{k)]e 


'.(lr\}*~ k ( Z+L ) 


-L 


- [e~ 2kL + l) J pdpJ,(kp)B.(pL) 


+ / 


ike~ lkL \p 2 dpJ x {kp)B (p (pL). 
0 


(36) 


i 

Now we multiply both sides of equation (35) by J poJ\{P m po)dpQ and get 

o 

y,(A) 


c » + "VO = 1 


*Pm A 


wi 


J * -A 


(37) 


The resulting expression involves a number of auxiliary matrices that can be evaluated numerically. 
These can be presented as follows: 


where 




J M r (o 

k 2 -Pl n 


(38) 


G 


(1)_ -2kL g ( 2)_ Ajk) L -2 kL\ Q (3) = kJ x {k)( 2kL\ 

’> J ’ U " i- 2 _v2V * i 


k ~Pn 


2 ' 
Yn 


G (4)_ Mk )_ 
~ , 2 2 

k -Yn 


\\ + e 


-2kL\ 


) 


^tS'l ^o(^;i)^l(^)~ [\[^Mjo{k) ( x -2kL\ 

" = + 1 1 
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(38a) 


G^ = -^{\ + e- 2kL ) , Gi 7) = ^§-( \ + e~ 2kL ) , c! 9 ) = ^(l-f 2W ) , 


** + P„ 




G (10)_ j i 4-^^ 

n 1 2 2 
A: — £ n 


(l + e - 2ii ),and = J 




i*(* 2 -& 2 ,) 


Now we define principal matrices: 


Yn 

'm,n . 3 


„ = ■ 


$m,n + 'o L { G ml + G mi] + tanh (7/r^) 


2fi Pn 


R q( 2 ) _XlLQ( 4 ) 
Hn^mn n w mn 

Pn j 


N m ,„ =(- 1 )" 


2 , 2 


c (5) + £a±5_f G (7)_ c «S)\_ , . g (6) 

w m n~ . 2 mn w m n "n m n 


_ ^0^(0 „ 

+ fi 0;(l) 

Pm 


,(■> + _!_ 

" 2A 2 


■ ^0) 


' 0 )_ 1 


nr J 


8 m 


wh 


(39a) 


(39b) 


■5Xf- 


/];(9) _/^(l0)« # an i 
w mn ^ mn °n iai1 


h(e„L)] , 


(39c) 


where 


*? 2/ 1 , Pji l ifej 


We end up now with the matrix equation 


[m]c = f+[n]d . 


(40) 


(41) 


Now consider the case of r = 1 , —L < z < L . We have from equation (34) 


-2*'flp(lzo) = J</*/i(*) n (*z 0 ) , 

o 


( 42 ) 
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where 


L 

n(£z 0 ) = - J dz[B r (\z)J { (k) + iB^\z)U[{k)]e k \ z - z ^ 

-L 


1 

+ e~ kL (e kz ° + e- kz °)\ pdpJ^kp^kpB^pL) - B._{pL )] . 

0 

In order to assure for a fast Fourier cosine series convergence, we need to extract the constant 
term from both sides: 


( 43 ) 


-2/[^,(ir 0 )-i*,(iz.)]= J</*y,(*)[n(fc 0 )-n(*L)] . 

0 

Then the cosine series expansion gives 

-D m -sk^- = J dz cos(p„,z)J rf«i(*:)[n(tzo)- n(*i.)] . 
Pm A m _ L Q 


(44) 


(45) 


After considerable algebra, this expression can be put into a matrix form. For convenience, we introduce 

a set of matrices, which can be numerically evaluated and are tjj) = J dkJ\ (k)r^ l \ where 

0 

m kjM ip ( 3 ) t 2 m) 

k 2 + pl’ *(* 2 +Pm) k 2 +p 2 m ' 

( 4 ) (' + ^ 2tL )Mk) ( 5 ) (l-e- 2it )./f« (6) k 2 (l + e- UL )j 2 (k) 

' = k(k 2 + P l) ’ ' = k 2 +P l = * 2 + p 2 ’ 


(7) k(l + e- 2kL )j’,(k) 

t n ^ 


k 2 +p 2 m 


, and 



(i -«*“)/,(*) 

[k 2+ Pm) 


(4.6a) 


and, similarly T^) n = J dkJ\[k + ) t^ , where 

0 
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r (') = 


1 

fl + e- 2W 

')JM 

(* 

2 +pl)( 

k 2 + Pn) 


, T 


(2) (i+‘- 2iL )tf!W 

(* 2 + Pm)(* 2 + pl) ' 


r (3) = 


k 2 

{l + e~ UL 

) Kio(KVM- 

■/, (A„)*/ 0 (t) 

( 

* 2 + Pm) 

I 

(* 2 + 
l n ) 

l'l(A„) 


4) (i +e - 2 “)P/,w (5) (l-KT™)* 2 /,**) 


_(«) (i+ g - 2tJ -)^ 2 J,W (7) (i- g - 2tt )* 3 j|(t) 

~(* 2 + />™)(* 2 - r ;;)’ ~ (* 2 + pj )(* 2 - r 2 ) ’ 


is) (i+«- 2 “y,(*)* 2 (9) (l-*- 2 * 1 )* 3 /,^) 

_ (* 2 + Pm)(* 2 - «n) ' a " '(* 2 + P 2 )(* 2 -e 2 ) 


We now introduce the following matrices: 


U - 2 L 
u mn 


^-'Ipj + q 2 | o) 

p„, -i 2 


W7 




- 2 */iPmH) 


m+/i+l 


Pm 


+ 2(-l) 


2 2 

m+n+1 Pn 


Pnjtt- — 

Pm 




PmPn 


_l£m!nl + T (4) j(7) t U 
2 ? ^ 1 mn ^ 1 mn ^ ia,u, V/ 


2 2 
7/i + Pm 


y 1 mn i n 1 mn 


(46b) 


(47a) 


(47b) 
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E«=*Or0)H) 


m + 1 


2 P 


,p) + ±. 

1 m • 


0 )' 

m 


m 1 m 


Pm 




m + 1 


4 M 3 )-l) 

2„ ,(5) + T w — — 2p t & 

Pm 


-Xw-ir!^ 




2? 


(l) 


2 2 
£/t Pm 


£ „tanh(£„L)r£>-r® 


(47c) 


Then the second matrix equation is 


[U]D + [V]C = E . 


(48) 


Combining this equation with equation (41), we obtain 


{[U]+[V][M] 1 [N]}D = E-[V][M] *F. (49) 

This matrix equation can be readily numerically inverted to yield the required values for vectors D and C. 

The required value C 0 can be obtained from the condition of zeroing the last C n coefficient. This condi- 
tion makes the solution nonsingular by selecting the proper value for the z -independent term. 
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6. NUMERICAL RESULTS 


In this section, we present the results of our numerical calculations based on the previous section. 
The aspect ratio is selected to be 1/2, which is close to our experimental value. The other relevant 
nondimensional parameter is \q\, which is proportional to the square root of the rotational speed. The 
dependence of the magnetic field on this parameter is needed for comparisons with the experiment. For 
illustration purposes, we present below a set of six graphs for three components of the field within the 
cylinder boundaries, each component having the real and the imaginary part (fig. 8). The graphs are 
given for q = V— 7 . The rank of the matrices entering equations (41) and (49) was selected to be 10. This 
number is rather small, but gives sufficient accuracy for this electromagnetic problem (10 -3 ). As expected, 
the B, component vanishes at the center of the cylinder. The B r value is largest at the middle of the rim, 
considering that we look only at the surface points. 



Figure 8. A set of six functions of r and z variables that represent the magnetic field 
in the spinning cylinder: (a) R e(#,.), (b) Im(£,.), (c) Re(£ y ), (d) Im(Ry), 

(e) Re(ZL) and (f) Im(£.). Re is the real part and Im is the imaginary part. 
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Our numerical analysis shows a high degree of diagonality of the matrix M. This was observed 
for practically all investigated L and q parameters. The nondiagonal elements are approximately five 
orders of magnitude smaller than the corresponding diagonal element. This feature can be exploited to 
obtain good analytical formulas, as the inversion of M can be accomplished analytically with high 
accuracy. Similarly, the governing matrix in equation (49) is highly diagonal as well (at least three orders 
of magnitude difference). Figures 9 and 10 represent the B r and fields taken at the middle of the rim 
at several different angular positions. As expected, the qualitative behavior is the same as in the case of 
an infinitely long cylinder. The difference is in the position of the maximum for angles in vicinity of 90°. 
Specifically, the peak shifts toward the higher q values. 




Figure 9. 5,./S 0 as a function of the skin 
depth parameter, q. Values are 
taken at the middle of the rim. 




Figure 10. B^/Bq as a function of the skin 

depth parameter, q. Values are taken 
at the middle of the rim. 


26 


7. EXPERIMENTAL RESULTS 


A brass (copper+zinc) wheel, with a 4-in diameter and 1 .4 in thick, was used for this investiga- 
tion. Ceramic ball bearings were used to safely spin the wheel up to 40,000 rpm. The wheel was also 
electrically insulated from the rest of the armature. For spinning, an air turbine was used. The rotational 
velocity was measured by two methods: (1) A simple strobe light and (2) a photodiode that detected 
laser beam pulses as they passed through a slot in the rotating shaft. The magnetic sensor is based on a 
giant magnetoresistivity, and consists of a balanced bridge circuitry. The position of the sensor was as 
close as possible to the rim of the wheel. The linear dimension of the sensor is =8 mm so that the offset 
from the surface is on the order of 1 5 percent. We did not use any goniometer system, so the accuracy of 
the angular position is not high, being estimated within a few degrees, with the main uncertainty being 
the direction of Earth’s magnetic field. We attempted to fit the experimental data with the presented 
theory by selecting the best value for the electrical conductivity of the wheel. The results of this proce- 
dure are displayed in figure 11, where the black dots represent experimental values. A slight misfit on 
the right shoulder can be due to slight angular misalignment from a 90° position. The obtained value for 
the resistivity is 43 nQm, which compares well with those listed in table 1 . We can conclude, based on 
these measurements, that the proposed theory satisfactorily explains our experiments. 



Figure 11. Experimental values (black circles) and the corresponding theoretical fit 
(hollow circles) assuming that the electrical conductivity is 43 nQm. 
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Table 1. Measured electrical resistivities for several brasses. 


Designation 

Composition 

Electrical 

Resistivity 

(n£im) 

C21000 

95 Cu - 5 Zn 

31 

C22000 

90 Cu - 10 Zn 

39.1 

C22600 

87.5 Cu - 12.5 Zn 

43 

C23000 

85 Cu - 15 Zn 

47 

C24000 

80 Cu - 20 Zn 

54 
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